Two-gap superconductivity in Bai_ x .K x .Fe 2 As 2 : A complementary study of the 
magnetic penetration depth by /iSR and ARPES 
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We investigate the magnetic penetration depth A in superconducting Bai_ :c K I ;Fe2As2 (T c — 
32 K) with muon-spin rotation (/iSR) and angle-resolved photoemission (ARPES). Using /iSR, we 
find the penetration-depth anisotropy 7a = A c /A a f, and the second-critical-field anisotropy 7ir c2 to 
show an opposite T-evolution below T c . This dichotomy resembles the situation in the two-gap 
superconductor MgB2. A two-gap scenario is also suggested by an inflection point in the in-plane 
penetration depth A a b around 7 K. The complementarity of fiSH and ARPES allows us to pinpoint 
the values of the two gaps and to arrive to a remarkable agreement between the two techniques 
concerning the full T-evolution of A a b. This provides further support for the described scenario and 
establishes ARPES as a tool to assess macroscopic properties of the superconducting condensate. 

PACS numbers: 76.75.-K, 74.70.-b, 74.25. Ha 



Much effort is devoted to the investigation of the man- 
ifestations and the mechanism of unconventional super- 
conductivity in the iron-arsenides, since many of their 
features clearly set them apart from other supercon- 
ductors. Ab-initio calculations, for instance, indicate 
that superconductivity originates in the d-orbitals of 
the Fe ion, which normally would be expected to be 
pair-breaking [l|, 0] • Several disconnected Fermi-surface 
sheets contribute to the superconductivity, as revealed 
by angle-resolved photoemission spectroscopy (ARPES) 
[1,0,11, El- Furthermore, indication for multi-gap super- 
conductivity was obtained in measurements of the first 
and second critical fields H c i and H c2 the magnetic 

penetration depth A [^, , as well as in point-contact 
Andreev reflection spectroscopy experiments 11 1 . 

Measurements of A provide a conclusive method to 
reveal multi-gap superconductivity, since the presence 
of gaps with different gap-to- T c ratios induces the ap- 
pearance of inflection points in A(T) [l2|, [II, [3, Ha |. 
Measuring A; (i = a, b, or c) along certain crystallo- 
graphic directions allows, in addition, the investigation 
of the penetration-depth anisotropy 7a . Within the Lon- 
don approximation, which implies A^ 2 oc n s /m* (n s is 
the carrier concentration), 7a is directly related to the 
anisotropy of the supercarrier mass m* via 7a = Aj/Aj = 

Jfn*j ffn* . As shown for the case of MgB2 [l6| , a differ- 
ent temperature evolution of 7a and the second critical 
field anisotropy 7_f/ c2 is also indicative of multi-gap su- 
perconductivity. 

Here we report a combined study of the penetration 
depth in a single crystal of Bai_ a; K a ;Fe2As2 (BKFA) by 
means of /iSR and ARPES. The sample was extensively 
characterized and several publications report its inves- 
tigation by ARPES, magnetic neutron scattering, /zSR 



and magnetic force microscopy [1, @, ^J. Resistivity 
and dc-susceptibility measurements demonstrate a sharp 
superconducting (SC) transition at T c = (32 ± 1) K, 
reproducible among different crystals from the same 
growth batch, and X-ra y p owder diffraction has estab- 
lished the phase purity 17 |. Most importantly for our 
study, the gap structure was investigated by ARPES 
Q. Furthermore, the occurrence of electronic phase sep- 
aration into antiferromagnetic (AF) and superconduct- 
ing/normal state regions on a lateral scale of several tens 
of nanometers was established [1, [13] • 

We begin by reporting on the /iSR measurements 
which were carried out at the 7rM3 beam line at the 
Paul Scherrer Institute (Villigen, Switzerland). The 
Bai_ x K x Fc2As2 single crystal with an approximate size 
of 5x10x0.06 mm 3 was mounted on a holder specially 
designed to perform fiSH measurements on thin single- 
crystalline samples. The transverse-field (TF) and the 
zero-field (ZF) /iSR experiments were performed at tem- 
peratures ranging from 1.5 to 200 K. In two sets of TF 
measurements the magnetic field (fioH = 10 mT > H c i) 
was applied in parallel and perpendicularly to the crys- 
tallographic c-axis, respectively, and always perpendicu- 
larly to the muon-spin polarization. The typical counting 
statistics were ~ 10 positron events for each particular 
data point. 

Experiments in transverse field allow to study the mag- 
netic ordering as well as to obtain the superfluid density 
response fljj ]. Muons stopping in magnetically ordered 
parts of the sample lose their polarization relatively fast, 
since the magnetic field on the muon stopping site be- 
comes a superposition of the external and the internal 
fields. The superconducting response is observed as an 
additional damping below T c because of the inhomoge- 



2 



neous field distribution of the external field penetrating 
the sample in form of vortices. Our ZF /iSR experiments 
reveal that the signal from the magnetically ordered parts 
vanishes within the first 0.3 fj,s. Bearing that in mind, 
in the whole temperature region the fit of TF data was 
restricted to times t > 0.3 fxs (see Ref. [H for details). 
For T < T c , the TF /j,SR data were analyzed using the 
following two-component form: 



A"(t) = A ie xp(-(^ 



Jt 2 /2) cos( 7 Bit + 4>) + 



A 2 exp(-cr^ m i 2 /2) cos{jB 2 t + 



(1) 



Here Ai and A 2 are the initial asymmetries of the first 
and the second component, 7/271 = 135.5 MHz/T is the 
muon gyromagnctic ratio, <p is the initial phase of the 
muon-spin ensemble, and the depolarization rates a sc and 
cr nm characterize the damping due to the superconduct- 
ing and the weak nuclear magnetic dipolar contributions, 
respectively. The second term on the right-hand side of 
Eq. ([I]) accounts for the parts of the sample remaining 
in the normal state [T3|. Each set of TF-/iSR data, con- 
sisting of measurements in the H \\ c and H T c config- 
uration, respectively, was fitted simultaneously with A±, 
cr nm , <fi, and B 2 as common parameters and cr sc , B\, and 
A 2 as individual parameters for each temperature point. 
The validity of our approach to fit some of the parame- 
ters globally was confirmed by examining the evolution of 
Ai, <Jnm and B 2 in the "free" fit. Above T c , the fit was 
simplified to the single Gaussian component only, with 
all parameters kept free. The results of the analysis are 
presented in Fig. [TJ 

The inset in Fig. [T] shows that the initial TF asym- 
metry A TF (t = 0) = A\ + A 2 (closed symbols) starts to 
decrease below T ~ 70 K, following the gradual enhance- 
ment of the magnetic fraction already investigated in this 
sample 17, 2(|. Here, we concentrate on the SC proper- 
ties and study in detail the temperature range belo w T f , 
which remained unexplored in the previous study [17j |. 
We note that A TF (closed symbols) and the asymmetry 
related to the superconducting fraction A\ (dashed lines) 
are almost the same for the H || c and H _L c sets of mea- 
surements. This is exactly what is expected, since these 
asymmetries must represent the corresponding volume 
fractions (magnetic or superconducting). 

The temperature evolution of the SC part of the muon- 
spin depolarization rate a sc is presented in the main 
panel of Fig. [T] It is worth to note that in a homogenous 
superconductor cr sc is expected to be proportional to the 
inverse squared magnetic penetration depth, er sc oc A~ 2 . 
In addition, in single-crystalline sample the magnetic 
field distribution in the superconductor in the mixed 
state is asymmetric and, therefore, can not be described 
by a single Gaussian line (see e.g. [HI, 14, [HI). The 
Bai_ 2 ,K a; Fe2As2 sample studied here is, on the contrary, 
highly inhomogeneous 17J, and the superconducting re- 
sponse, at our level of statistics, is well described by a 
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FIG. 1: (Color online) Temperature evolution of a sc measured 
after field cooling in noH = 10 mT applied in parallel (red cir- 
cles) and perpendicularly (blue circles) to the crystallographic 
c-axis. The inset shows the temperature evolution of the ini- 
tial TF asymmetry A TF (t = 0) = Ai + A2 (closed symbols) 
and the superconducting asymmetry A\ (dashed lines) nor- 
malized to ,4 TF at T = 100 K. The colored areas represent 
volume fractions. Note the logarithmic T-scale. 



single line of Gaussian shape. We believe, however, that 
for Bai_ a: K 2 ,Fe2As2 er sc is still a good measure of A. In- 
deed, cr sc at H || c extrapolated to T — * results in 
cr sc (0) ~ 1.2 /us -1 , which follows reasonably well the Ue- 
mura relation established recently for various families of 
Fe-based superconductors 21 1. We conjecture that 
the antiferromagnetic islands act as preformed pinning 
centers for vortices, thus precluding the formation of an 
ordered vortex lattice [l7| , while the screening current at 
this relatively low field (/j,qH = 10 mT) still flows at a 
distance A from the vortex core. 

Within the London model, the magnetic penetration 
depth of the isotropic extreme type-II superconductor 



(A ^> £, £ is the coherence length) is determined by A oc 
n s /m*. For an anisotropic superconductor, the magnetic 
penetration depth is also anisotropic and is determined 
by an effective mass tensor [22]: 



Mi 



<f = Mj 



M k 



where M, 



and 



(2) 



is the mass of 



n*/f/m*-m* • „„ fc 
the carriers flowing along the i-th principal axis. The 
effective penetration depth for the magnetic field applied 
along the i-th principal axis of the effective mass tensor 
is then given by [12| : 



A eff — 



A, A 



(3) 



For convenience, we drop the index "sc" in the "super- 
conducting" Gaussian relaxation rate er sc . Equation Q 
implies that by applying the magnetic field along the 
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FIG. 2: (Color online) Temperature evolution of the magnetic 
penetration depth anisotropy 7a = A c /A a i, = cr" c /er Xc . The 
line is a guide to the eye. In the inset we compare 7a with the 
i/ C 2-anisotropy "fH c2 obtained for Bai_ a; K a ;Fe2As2 by Yuan 
et al. [23[ and Altarawneh et al. [24|, albeit in samples with 
somewhat different T c . For better comparison, T is divided by 
the respective T c and the values of 7a and 7« c2 are normalized 
to 1 for the data point closest to T c . 



MeFeAsOi_ x F x (Me=Nd and Sm) ^ and the double- 
layer BaFc 2 _ x Co x As2 [13]. The authors of Ref. [H also 
explain the observed behavior by the presence of multiple 
gaps opening on various bands at the Fermi level. 

An additional confirmation for the multi-gap behavior 
comes from the analysis of the temperature dependence 
of the in-plane magnetic penetration depth. Fig. [3] shows 
A~ b 2 (T) obtained from the measured a^ c (T) by using the 
relation cr^s -1 ) = 0.1067 • A" 2 (/im -2 ) (28|. The experi- 
mental A~ b 2 (T) data were analyzed within the framework 
of the phcnomcnological a-modcl by assuming two inde- 
pendent contributions to the total A~ h 2 [l^ |: 

KbiT) = A^ a (0)( W .{l-D[Ai(T),T]} 

+(l-w)-{l-D[A 2 (T),T]}), (6) 

where D(A,T) = f_ + ~ Re-^de [13, /(e) 

is the Fermi function, and cu is the contribution of the 
bigger gap to A~ fc 2 . The temperature dependence of the 
superconducting gap is assumed to be [12[ 

Ai, 2 (T) = A li2 (0)tanh{1.82[1.018(T c /T- l)] - 51 }, (7) 



crystallographic a, b, and c directions, one measures 
<j" a oc l/At,A c , er" b oc l/A a A c and er" c oc l/A a Af,, respec- 
tively. By neglecting the difference between A a and A^ 
the penetration depth anisotropy can be obtained as: 



7a = A c /A o6 = (T l|c /cr J 



(4) 



The temperature evolution of 7a is presented in Fig. [5] 
The inset shows the second critical field anisotropy 
Jh c2 = H^f/H^ obtained for similar Bai_ a; K a; Fe2As2 
samples in resistivity [23| and radio frequency penetra- 
tion depth measurements (24j. Within the phenomeno- 
logical Ginzburg-Landau theory, in a single-gap super- 
conductor both anisotropies must be equal |25| |: 



7a = 




= 7ffc2 



Cab 



It is natural to expect that the values of the same 
quantities measured by various techniques should be the 
same. While a deviation at a particular temperature 
might be explained by a slight variation of the properties 
among samples used in the different experiments, this 
cannot account for the opposite temperature evolution 
of 7a and jh c2 shown in the inset of Fig. [2 Hence, in 
Bai_ a; K a ;Fc2As2 7a and 7_y c2 are not the same and Eq. (|5|) 
is violated. This resembles the situation in the two-gap 
superconductor MgB 2 , albeit with reversed trends for 7a 
and 7/f c2 : In MgB 2 , 7a was found to decrease with de- 
creasing temperature from about 2 to 1.1, while 7// c2 
increases from ~2 at T c to 6 at low temperatures [3]. It 
is worth noting that recently the presence of two distinct 
anisotropics 7a and Jh o2 was reported for the single-layer 



in agreement with A(T) measured by ARPES (see 
Fig. Eh). The solid line in Fig. [3^ represents the result of a 
fit of Eq. ((6]) to the experimental data with A a b(0), w, Ai , 
and A2 as free parameters. The fit yields: Ai =9.1 meV, 
A 2 = 1.5 meV, lu = 0.5, and X ab {0) = 320 nm 

The penetration depth A a b(T) can also be calculated 
from the electronic band dispersion and the momentum- 
resolved SC gap [3l[ which were determined by ARPES 
on BKFA single crystals from the same growth batch 
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FIG. 3: (Color online) Temperature evolution of the inverse 
squared in-plane magnetic penetration depth obtained 
from the measured cr" c presented in Fig. [l]by using the rela- 



tion o-(^s _1 ) = 0.1067A" 2 (/im" 



The solid line repre- 



sents the result of a fit of Eq. (0 to the a-model, the dashed 
line represents a calculation of A~ 6 2 from the electronic struc- 
ture revealed by ARPES [2^] with one fitting parameter, A2 . 
Inset: contributions of different Fermi surface sheets to A~ b . 
(b) Fermi surface of BKFA. (c) Temperature dependence of 
the SC gap, extracted from ARPES spectra [f|. 



4 





/ibK 


Axvr £jo 


Aab(O) (nm) 


320 


270 


LU 


0.51 


0.55 


Ai (meV) 


9.1 


9.1 


A 2 (meV) 


1.5 
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TABLE I: Parameters as extracted from the fit of /iSR data 
and calculated from ARPES data. 



@, H El, H3] ■ The Fermi surface consists of four differ- 
ent sheets — an inner T-barrel, an X-pocket and blade- 
shaped pockets with a large isotropic gap Ai, and an 
outer T-barrel with a small gap A2 |fj, |32j (Fig. [3b) ■ The 
formula relating A to the electronic structure reads 

Kb( T ) =h{l- D[Ai(T), T]} + I 2 {1- D[A 2 (T),T]} , 

(8) 

where I12 are integrals over the Fermi-surface contours 



h = 



«p(k)dfc, h = 



«F(k)dfc, 



2irsoc 2 hL c J ' 2ireoc 2 hL l 

outer r, inner T 

blades 
X-pockct 

(9) 

£0, h, e, c are physical constants, L c is the c-axis lat- 
tice parameter, and v-p is the Fermi velocity. Further 
details of the electronic structure and the calculation are 
given in ref. [2^. Eq. ([8]) is equivalent to Eq. (® with 



Kb (°) = h+h and a; 



V6 \"> ~ ^ ' ^ ^ - T[+h - Usin S Ec L ©' we Cal " 

culate A~ b 2 (0) and 10, while Ai is known from the ARPES 
measurements [H|, @ . A comparison of these parameters 
determined by the two different methods is shown in Ta- 
ble I. Taking into account the complementarity of the 
methods, the agreement is remarkable and strengthens 
the validity of the obtained results. The discrepancy in 
A a f,(0) is well within the range of values obtained in dif- 
ferent //SR experiments (icl | . 

Now we can assess the remaining parameter A2 by 
fitting Eq. ([5]) to the measured X a b(T) normalized to A(0). 
fn Fig.[3]the normalization is realized by scaling the A~ SR - 
and A^ PES -axes accordingly. We obtain A 2 = 1.1 meV, 
again in good agreement with the /iSR result. 

In summary, from /zSR measurements on a single- 
crystalline sample of BKFA (T c ~ 32 K) we have deter- 
mined the anisotropy of the magnetic-field penetration 
depth 7a = X c /Xab- The penetration depth anisotropy 
increases with decreasing T from 7 a ~ 1.1 at T ~ T c 
to 7a ~ 1.9 at T ~ 1.7 K, while the T-evolution of the 
H C 2 anisotropy Jh o2 shows an opposite trend [2^, 24 1. 
This resembles very much the situation in double-gap 
MgB 2 where both anisotropies are equal at T c , but evolve 
oppositely with T. The notion of two SC gaps is sup- 
ported by the observation of an inflection point in A Q b 
at ~ 7 K. From a fit of A~ 6 2 to the phenomenological 
a-model we obtain gap values of Ai = 9.1 meV and 
A2 = 1.5 meV. A comparison of A~ h 2 (T) measured by 



/iSR with the one calculated from ARPES data shows a 
remarkable agreement between these two complementary 
approaches, lending further support to our conclusions 
and establishing ARPES as a tool to estimate A a b. 

The fj.SK work was performed at the Swiss Muon 
Source (S^S), Paul Scherrer Institute (PSI, Switzerland). 
ARPES results were obtained at BESSY. 
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